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ABSTRACT 

Fourier decomposition is a well established technique used in stellar pulsation. 
However the quality of reconstructed light curves using this method is reduced 
when the observed data have uneven phase coverage. We use simulated annealing 
techniques together with Fourier decomposition to improve the quality of the 
Fourier decomposition for many OGLE LMC fundamental mode Cepheids. This 
method restricts the range that Fourier amplitudes can take. The ranges are 
specified by well sampled Cepheids in the Galaxy and Magellanic Clouds. 

We also apply this method to reconstructing Cepheid light curves observed 
by the HST. These typically consist of 12 V band and 4 I band points. We 
employ a direct Fourier fit to the 12 V band points using the simulated annealing 
method mentioned above and explicitly derive and use Fourier interrelations to 
reconstruct the I band light curve. We discuss advantages and drawbacks of 
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this method when apphed to HST Cepheid data over existing template methods. 

Application of this method to reconstruct the light curves of Cepheids observed 
in NGC 4258 shows that the derived Cepheid distance (/xq = 29.38 ±0.06, random 
error) is consistent with its geometrical distance (/Xq — 29.28 ±0.09) derived from 
observations of water maser. 

Subject headings: Cepheids — methods:data analysis — distance scale 



The study of Cepheid light curves has two major applications: (a) distance determina- 
tions to nearby galaxies via the well calibrated period-luminosity (PL) relations (for example, 
see Madore & Freedman (1985); Feast & Walker (1987); Madore & Freedman (1991); Saha 
et al. (2000c)); (b) the understanding of the pulsational behavior of Cepheids by comparing 
the observed light curves to the theoretical counterparts (Simon & Davis 1983; Buchler et 
al. 1990). Therefore, reconstructing the Cepheid hght curves from the observed photometric 
data is important regarding these two aspects. 

Because Cepheid light curves are periodic, the data points from well observed Cepheids 
can be described by Fourier expansion. This technique was first introduced by Schaltenbrand 
& Tammann in 1971, and developed by Simon & Lee (1981) to study the structural properties 
of Cepheid light curves. In general, the n*'*-order Fourier expansion has the following form: 



where u; = 27r/P is the frequency and P is the period of Cepheid. The Aq term is the 
mean value of the light curve; the Aj and (f)j are the Fourier amplitudes and phases for j^^ 
order, respectively. In this paper we develop the use of simulated annealing techniques to 
improve the quahty of light curve reconstruction in nearby galaxies such as the LMC. We 
also apply these techniques in the light curve reconstruction of Cepheids observed by Hubble 
Space Telescope (HST) which is necessary for distance determinations. In particular, we 
use Fourier expansion, as given by equation (1) but with modifications, to reconstruct the 
V band light curves and Fourier interrelations to reconstruct the I band light curves. These 
Fourier techniques have been established in stellar pulsation studies. We also present some 
potential advantages and drawbacks of the application of our method to this problem. 



1. 
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Section 2 & 4 will describe the Fourier expansion and interrelations in details, respec- 
tively. Section 3 briefly describe the application of these Fourier techniques to HST data. 
The error analysis of the light curve reconstruction procedures will be presented in Section 
5, following by the discussion and conclusion in Section 6. In appendix, we show another 
Fourier method, the Fourier intrarelations, that can be applied to reconstruct the Cepheid 
light curves. However, this method is less preferable than Fourier interrelations. 



The form of Fourier expansion is given by equation (1), which has been extensively 
used in studying the light curves of pulsating stars (see Antonello et al. (1987) and reference 
therein). Some applications of equation (1) in the study of the Ccphcids can be found, for 
example, in Simon & Lee (1981); Simon & Davis (1983); Simon & Moffett (1985); Antonello 
& Poretti (1986); Andreasen & Petersen (1987); Buchler & Moskahk (1994); Poretti (1994); 
Hintz & Joner (1997). Obviously, equation (1) has 2n + l unknown parameters which require 
at least the same number of data points to solve for these parameters (in this paper we assume 
all of the periods for the Cepheids are given in the literature because they can be determined 
with other methods). Most of the observations on Galactic (eg., Schaltenbrand & Tammann 
(1971); Simon & Lee (1981); Moffett & Barnes (1985); Ferro et al. (1998)) and Magellanic 
Clouds Cepheids (eg., Andreasen & Petersen (1987); Moffett ct al. (1998); Udalski et al. 
(1999a,b)) include large number of epochs that permit the higher order of Fourier expansion 
up to n = 8. In contrast, there are typically about 12 V band epochs in the HST observations 
of extra-galactic Cepheids (see Section 3), which permit a 4*'* order Fourier expansion. The 
4*'* order fit is a compromise between the goodness of fit and available data: higher order fits 
may not result in a better fit and lower order fits may not capture the structural properties of 
the light curves. Since there are only 4 I band data points in HST observations (see Section 
3), the 4*'* order Fourier expansion cannot be applied. 

Since the period is known from the literature, we can fold the time observation into 
phase as: ^{t) — {t — to)/P — int[{t — to)/P], with to being a common starting epoch. The 
value of $ is from zero to one, corresponding to a full cycle of pulsation. Hence equation (1) 
can be written as (Schaltenbrand & Tammann 1971): 



In general, one approach to solve the unknown parameters in equation (2) is by using 



2. Fourier Expansion 
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the linear least-square method together with the available data. However, we find that in 
some cases the solution does not provide a good fit to the data and produces some numerical 
bumps in the light curves. Figure 1 illustrates this problem. The solid curve is the resulting 
light curve from the least-square solution to a Cepheid in NGC 1326A. The light curve 
does not fit to the data well and exhibits numerical bumps around phase ~ 0.2 and ~ 0.5. 
These numerical bumps are not physically associated with the true light curves, but result 
Irom the sparsely sampled data points. Instead of solving for the unknown parameters, we 
specify the range of each parameter in equation (2) and fit the data to obtain the best-fit 
values for each parameter. All the parameters are fit simultaneously with the simulated 
annealing method (for example, see Press et al. (1992)) to minimize the corresponding 
value: = S[(^o6s — '>TT'fit)/c'obs]'^- The simulated annealing method will search the best 
fit values for a group of parameters within a reasonable amount of computational time. 
The dashed curve in Figure 1 shows that by fitting the parameters in equation (2) to the 
observed data with an appropriate choice for the ranges (given in Section 2.2 for Fourier 
amplitudes, Aj), the reconstructed light curve has been greatly improved. In the following 
subsections, we demonstrate how the choice of the ranges for Fourier parameters will affect 
the reconstructed light curves. 



2.1. The Ranges of the Fourier Pcirameters 

In order to apply the simulated annealing method to fit the available data points, certain 
ranges of the fitted parameters have to be imposed. The range for the mean magnitude {Aq) 
term in equation (2) is trivial (say, ±5mag from the mean magnitude of the data points), 
and the range for the phases {(f)j) is between and 27r. For the Fourier amphtudes, the 
following ranges are initially used to fit the Fourier expansion: 

Al : 0.0 - 0.7; A2 : 0.0 - 0.5; A3 : 0.0 - 0.5; A4 : 0.0 - 0.5 

These initial ranges are selected based on the guess that some Cepheids may have large 
Fourier amplitudes, but Simon & Lee (1981) show the Fourier amplitudes are generally 
smaller than 0.4. Therefore we pick the above ranges as a starting point. By fitting the data 
with the ranges given above, most of the light curves can be reconstructed fairly well with 
good agreement to the data points. Figure 2(a) shows an example of the reconstructed light 
curve for a Cepheid in the LMC with this initial choice for the ranges of Fourier amphtudes. 
However, these initial ranges also produce some poorly reconstructed light curves, due to the 
bad phased coverage of the data points, as illustrated in Figure 2(b). Comparing the two 
cases in Figure 2 (which have similar periods) , it is clear that some of the Fourier amplitudes 



- 5 - 



for poorly reconstructed light curves are relatively large, indicating that the choice of the 
initial ranges is inappropriate. 

The quality of the reconstructed light curve can be improved by further narrowing the 
ranges of the Fourier amplitudes when fitting equation (2) to the data points. In order to 
determine the proper ranges for the Fourier amphtudes, we use the OGLE LMC data (Optical 
Gravitational Lensing Experiment, Udalski ct al. (1999a)) which is obtained from the OGLE 
home-page^. There are 771 fundamental mode Cepheids in the data, based on the judgment 
made by OGLE team (Udalski et al. 1999a). The number of epochs in V and 1 band are 
12 — 50 and > 80, respectively. The large numbers of epochs are sufficient to construct the 
light curves for individual Cepheids by using equation (2). However, the 10 Cepheids which 
do not have V band data and the Cepheids which are not fundamental mode Cepheids were 
removed from the final sample. Because the Fourier amplitudes are independent of distance, 
we also include a "calibrating set" of Cepheids (Hendry et al. 1999), which comprise mostly 
of Galactic and some LMC/SMC Cepheids. In addition, most of the Cepheids observed 
by OGLE team have periods shorter than 10 days, with peak around 0.5 < log(P) < 0.7, 
and terminated at log(P) = 1.5. In contrast, the "calibrating set" Cepheids occupy a wider 
range of period, from log(P) 0.2 to log(P) ~ 2.0. Therefore, including the "calibrating set" 
Cepheids can be used to determine the ranges of Fourier amplitudes beyond log(P) = 1.5. 
There are about 118 Cepheids in this sample. Since I band data for OGLE LMC has a large 
number of epochs, we mainly impose ranges for the Fourier amplitudes in V band. 

The V Band Fourier Amplitudes: After constructing the light curves of OGLE 
LMC Cepheids with initial (and relatively large) Fourier amplitudes, we divided the V band 
light curves into two groups by visual inspection: those with well constructed light curves 
(as in Figure 2(a)) and those with poorly constructed light curves (as in Figure 2(b)). The 
distributions of the Fourier amphtudes as a function of log(P) are plotted in Figure 3. 
This shows the positions of Fourier amplitudes for well constructed light curves (in crosses) 
and poorly constructed light curves (in triangles), as a function of period. As can be seen 
from the figure, the distributions of the Fourier amplitudes for well constructed light curves 
occupy a smaller range than the poorly constructed light curves, which have relatively larger 
amplitudes especially for higher order harmonics. The V band light curves are further 
trimmed to contain only the really well constructed light curves, exemplified by Figure 2(a). 
We then re-plot the distributions in Figure 4, which shows evidence of well defined ranges 
for the V band Fourier amphtudes. 

The I Band Fourier Amplitudes: Since the number of epochs for most of the OGLE 
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LMC Cepheids in the I band is more than 80, the constructed hght curve are expected to be 
well behaved and the Fourier amplitudes are assumed to fall within the appropriate ranges. 
The plots of the Fourier amplitudes are shown in Figure 5, with some suspected outliers 
marked in triangles and labelled. We also plot out some of the Fourier light curves for these 
outliers in Figure 6. As seen from the figure, some of them have acceptable light curves, 
with slightly larger amplitudes when compared to the Fourier light curves from Cepheids 
with similar periods. The exception is C686, which exhibits a tip around phase 0.1 from the 
original data. The simple 4*'* order Fourier expansion cannot reproduce this tip and hence 
the constructed hght curves do not show good agreement to the data. Nevertheless, the 
difference in amplitudes for the outliers and the majority of the data points are very small, 
because the I band amplitudes for Cepheids are generally smaller (Freedman 1988). 



2.2. The Adopted Ranges for Fourier Amplitudes 

From Figure 3-5, it is clear that the Fourier amplitudes occupy certain ranges in the 
amplitude vs. log(P) plots. Therefore we can determine the appropriate ranges of Fourier 
amplitudes for a given period from these figures. The adopted ranges for the Fourier am- 
plitudes in V and I band are presented in Table 1. Due to the fact that the OGLE LMC 
data only contains Cepheids with periods less than 32 days (log(P) < 1.5), the ranges for 
Cepheid with log(P) > 1.5 arc determined from the "calibrating set", while the ranges for 
shorter period Cepheids are determined from both OGLE LMC and "calibrating set" data. 
In addition, the ranges of the Fourier amplitudes are set to be started from zero to account 
for possible low amplitude Cepheids. Note that the upper limit of the ranges given in Table 1 
is an approximation, since there is no exact upper limit for a given period. Experience shows 
that sometimes a slightly larger range of the Fourier amphtudes than the one given in Table 
1 can reconstruct the light curve better. An iterative process can be made, if necessary, to 
find the most suitable upper limits to reconstruct satisfactory light curves. The reduction 
in ranges can improve the quality of the Fourier fit, as some of the numerical bumps arc 
removed in the reconstructed light curve, provided that the problem of bad phase coverage 
is not too severe (see details in Section 2.3). 

The distributions of the Fourier amplitudes for OGLE LMC Cepheids and the "cali- 
brating set" (or mostly Galactic) Cepheids appear to coincide, as seen in Figure 3-5. This 
may imply that the ranges of the Fourier amplitudes depend weakly on mctallicity. How- 
ever, this conclusion is only based on the analysis of two galaxies, and may not reflect the 
assumption that metallicity can affect the distribution of Fourier amplitudes, van Gen- 
deren (1978) showed that the upper hmits of B band amplitudes are different for Cepheids 
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in Galaxy/M31, LMC and SMC. In addition, Paczynski & Pindor (2000) showed that the 
OGLE LMC Cepheids, in the period range of 1.1 < log(P) < 1.4, have larger amphtudes 
than OGLE SMC Cepheids with same period ranges. Because the V band amphtudes can 
be scaled from B and I band amplitudes (Freedman 1988), the different V band upper limits 
for Cepheids in different galaxies could also exist. Nevertheless, the ranges of Fourier ampli- 
tudes given in Table 1 are assumed to cover the Fourier amplitudes for different metallicity 
environments, because these ranges go from zero to the upper limit that is slightly larger 
than the ranges defined by OGLE LMC and "cahbrating set" Cepheids. The detailed analy- 
sis of the relationship between the Fourier amplitudes and the metallicity environments will 
be presented in future work. 



2.3. Results and Problems of the Fits 

Using the techniques of 4*'* order Fourier expansion with the appropriate ranges for the 
Fourier amplitudes, some of the poorly reconstructed light curves can be improved and the 
results show good agreements to observational data. Figure 7 shows some examples of the 
improved reconstructed hght curves when fitting equation (2) to the OGLE LMC V band 
data by narrowing the Fourier amplitudes as given in Table 1. As can be seen from the 
figure, some of the obvious bumps (or dips) in the original reconstructed light curves can be 
removed with this simple technique, resulting in smoother light curves and better fits to the 
data. The estimation of the mean magnitudes would be closer to the true mean without the 
influence of the numerical bumps/dips. 

In certain extreme cases that cither the data points are heavily clustered in certain 
phases or there is a large gap between two adjacent data points, the Fourier expansion fails 
to construct a satisfactory light curve, regardless of the choice of the ranges for Fourier 
amphtudes. Figure 8 shows some examples of this case. By using either the initial ranges 
(dotted curves) or the narrower ranges (solid curves) to fit the Fourier amplitudes with 
observed data points, the reconstructed light curves are still not satisfactory due to the 
phase clustering of the original data. When using this technique to estimate means for the 
PL relations, these Cepheids should be treated with caution or rejected. In fact, the rejection 
of Cepheids with bad phase coverage has been practiced in the past (for example, see Simon 
& Lee (1981); Antonello & Poretti (1986)). 
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3. Application to HST Data 

Since Ccpheids arc primary distance indicators, and the Ccpheid PL relations can be 
used to calibrate secondary distance indicators, such as the Tully-Fisher relation. Type 
la supernovae, fundamental plane and surface brightness fluctuations. This builds up the 
distance ladder and indicates the distance to more remote galaxies. Together with the 
measurements of the recession velocities of these more distant galaxies, the Hubble constant 
can be determined: this was the main goal of Hubble Hq Key Project (Preedman et al. 
(2001), hereafter HOKP). 

The conventional way of applying the Cepheid PL relations is by measuring the mean 
magnitude of the Cepheid over the pulsational cycle for a given period. In recent years, a 
great deal of effort has been made to discover extra-galactic Cepheids, and determine the 

distance, by using the HST (as in HOKP and Saha et al. (1999, 2001a,b)). In most of the 
HST observations, in order to optimize the phase coverage of Cepheids within the HST 
observational windows, the number of V band observations is chosen to be 12 (Preedman 
et al. 1994; Kennicutt et al. 1995). To correct for the extinction and reddening, I band 
observations are included (see, for example, HOKP and Kennicutt et al. (1995)). However, 
only 4 (or 5) I band observations are needed, because the I band amphtudes are smaller 
than in the V band (Preedman 1988; Preedman et al. 1994; Kennicutt et al. 1995). Existing 
techniques to estimate the mean from the sparsely sampled HST data include: (1) taking 
the phased weighted intensity mean (Saha & Hoessel 1990) in V band and using an empirical 
relation developed by Preedman (1988) in I band (for example, see Silbermann et al. (1996)); 
(2) adopting template fitting procedures (Stetson 1996; Tanvir et al. 1999). 

In the application of the methods developed in the previous sections to this problem, we 
use we use Pourier expansion, as given by equation (1) but with modifications, to reconstruct 
the V band light curves and Pourier interrelations to reconstruct the I band light curves. 
Some results of applying these methods to determine the Cepheid distance in nearby galaxies 
are presented in Section 6.4. 



4. Fourier Interrelations 

Because the I band observation of HST typically consist of 4 epochs, equation (2) cannot 
be applied to reconstruct the light curves. However, statistical relationships between the j*'* 
order Pourier coefficients for V and I band have been introduced by Hendry et al. (1999) 
and by Kanbur & Nikolaev (2001), known as the Pourier interrelations. We explicitly derive 
and present them here because of their applicability in stellar pulsation studies. These allow 
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the reconstruction of I band light curves up to 4*'* order. There also exist correlations of 
the first order Fourier amplitude (Ai) to subsequent higher order Fourier amplitudes in the 
same bands (for example, see Antonello et al. (1987)), known as the Fourier intrarelations 
(presented in Appendix), this approach is less preferable than Fourier interrelations. The 
Fourier interrelations are the linear relation of j*^ order Fourier amplitude and phases from 
V band to 1 band, and have the following form: 



The coefficients of the interrelations are determined from the "calibrating set" , by min- 
imizing the value of the fit. Fourier interrelations for fundamental mode Cepheids are 
presented in Table 2. The table also shows the rms of the straight-line fit to the data. The 
fits were derived taking full account of errors in both variables with the standard model 
fitting procedures (for example, see Press et al. (1992) Chapter 15). The Fourier interrela- 
tions for equation (3) & (4) are clearly shown in Figure 9 & 10, respectively. The standard 
errors of the Fourier coefficients, derived from inverting the Hessian matrix, are also shown 
in the figures. Despite the appearance of strong correlations in the figures, the relations have 
large P^r degree of freedom, i.e. very small statistical significance. Nevertheless, using 
interrelation to recover sparsely sampled light curves works quite well, as shown in Section 
6. For the Galactic data in the "calibrating set", we display bump and non-bump Cepheids 
with solid and open squares respectively (due to the P2/P0 — 0.5 resonance around 10 days 
on Hertzsprung progression, the periods of bump Cepheids normally lie in between 8 and 
14 days.). For the LMC/SMC data, we do not differentiate between bump and non-bump 
Cepheids. 

Besides the Fourier interrelations among the "calibrating set" Cepheids, similar hnear 
relations also exist for OGLE LMC and SMC Cepheids. The Cepheid data of OGLE LMC 
are as in Section 2, and the data of OGLE SMC (Udalski ct al. 1999b) arc downloaded from 
OGLE home-page. The number of the fundamental mode Cepheids, as judged by Udalski 
et al. (1999b), observed in SMC is ~ 1300, with most of them in short period ranges. Table 
3 presents the best-fit results of Fourier interrelations to the OGLE LMC data, and the 
plots of the Fourier interrelations are given in Figure 11 & 12 for Fourier amplitudes and 
phases, respectively. The prehminary results of the OGLE SMC Fourier interrelations are 
presented in Table 4, with the plots of Fourier amplitudes and phases in Figure 13 & 14, 
respectively. Again a clear linear relation exists with increasing scatter as we go to higher 
order parameters. There is little difference in the slopes of the best fit lines to that found for 
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"calibrating set" Cepheids. Hence, the Fourier interrelations given in Table 2 will be applied 
to reconstruct the I band light curves. 

The procedures for reconstructing the I band hght curves are similar to V band hght 
curves. Instead of fitting all the parameters (the mean magnitude and the Fourier parame- 
ters) , we use the Fourier interrelations to obtain the I band Fourier parameters from the V 
band fit. This gives the shapes of the I band light curves, and the mean I band magnitudes 
are found by using the observed I band data to minimize the corresponding -y^ values. 

There is little difference between the Fourier interrelations for Cepheids in "calibrat- 
ing set" and in OGLE LMC or SMC, as shown in Table 2, 3 & 4, respectively. The good 
agreement between them indicates that the relative changes in Fourier parameters are al- 
most unaffected by metallicity. This does not mean, however, that the light curve shape is 
independent of metallicity. Rather, it means that change in metallicity affects the Fourier 
parameters in different bands in the same way. Therefore, given a well-sampled light curve 
in V band, one can estimate Fourier parameters in I band reasonably well regardless of the 
metallicity of the parent galaxy. This result is useful in reconstructing the light curves of 
extra-galactic Cepheids in a broad range of metallicity environments. 

The empirical relation that the earlier papers of HOKP used to obtain the I band mean 
magnitude is based on the observational results of Freedman (1988), as the amplitude ratio 
of I and V band is about 0.5. However, Tanvir (1997) recommended use of an amplitude ratio 
of 0.6, because this ratio can improve the estimation of I band mean magnitude. We note 
that the slopes in the Fourier interrelations {J5 term in Table 2 & 3) are approximately equal 
to 0.6, which is close to the Tanvir's value. This is the main reason why Fourier interrelation 
exists and work well in reconstructing the light curves (Section 6). 

5. Error Analysis for Light Curve Reconstruction 

The Fourier techniques described in this paper are now applied to reconstruct the 
Cepheid hght curves in nearby galaxies observed by HST, particularly for HOKP galaxies. 
Perhaps the most important question regarding the application of our work to this problem 
is the following: does a fourth order fit with 9 free parameters to 12 V band data points 
represent anything physically? In order to estimate the errors associated with the light curve 
reconstruction procedures presented in Section 2 & 4, we performed an error analysis based 
on Monte-Carlo simulations. 

Due to the large number of Cepheids, as well as the large number of observations in V 
and I bands, the simulations are performed primarily with OGLE LMC Cepheids that have 
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good light curves (as in Figure 2(a)) and represented by crosses in Figure 3. These OGLE 
LMC Cepheid hght curves, constructed by using all the available data points in both bands, 
are referred as the original light curves. We assume that when we fit a 4^^ order Fourier 
expansion to these data, the resulting Fourier amplitudes are very close to their true values. 

In order to mimic the published HST photometric data, we performed three simulations 
for the error analysis. The procedures for each simulation are as follows: 

• Simulation 1: First a Cepheid is picked randomly (with replacement) from the data 
set, then 12 points in the V band and 4 points in the 1 band were randomly selected 
without replacement. In addition to the OGLE LMC photometric errors, we add 
Gaussian noise of Unoise = 0-05mag and anoise = 0-lOmag to these randomly selected 
data. 

• Simulation 2: Same as in Simulation 1, but with larger Gaussian noise of /inoise — 
O.lBmag and keep anoise — O-Wmag. 

• Simulation 3: We pick one Cepheid with large number of epochs {N^ — 34, Nj — 
188), then 12 points in the V band and 4 points in the I band were randomly selected 
without replacement form this Cepheid. The additional Gaussian noise is same as in 
Simulation 1. 

Then the data from each simulation arc used to reconstruct the light curves with the 
procedures described in Section 2 & 4, i.e. 4*'* order Fourier fit in the V band and in- 
terrelations in I band. The randomly selected 12 V points could of course be uniformly 
distributed or concentrated around one particular phase point. Mean magnitudes and the 
Fourier amplitudes obtained from our reconstruction procedures for these simulated data 
are then compared to the mean magnitudes and Fourier amplitudes from the original light 
curves. This is repeated a large number of times to build up an error distribution. 

We emphasize that we are simulating the photometric data published in the literature. 
It has been suggested by an anonymous referee that even this is not a real simulation of 
HST Cepheid data since we neglect the possibility of events such as "warm pixels" or cosmic 
rays. However we contend that such data points will either be rejected by the photometric 
reduction package used and, if not, the point sources responsible for these should not be 
used anyway. 

After running N — 1000 trials in each simulations, the error histograms for the mean 
magnitudes, as well as the Fourier amplitudes, in both V and I bands are constructed. Gaus- 
sian distributions with parameters of /i and a are then fitted to these error histograms, where 



- 12 - 



represent the mean offsets between the simulated data and original data, and a represent 
the errors in either the means or the Fourier amplitudes. The results of the Gaussian fits 
from each simulations for V and I bands are presented in Table 5 & 6, respectively. The error 
histograms for the mean values, resulted from Simulation 1, are presented in Figure 15, with 
abscissa to be the difference between the simulated means and original means. Similarly, the 
errors histogram for Fourier amplitudes in V and I bands are presented in Figure 16 & 17, 
respectively. The parameters from Gaussian fits are hsted in the upper-left corners of these 
figures. We did not include the error histograms for Simulation 2 & 3 in this paper, because 
they all look similar to the error histograms of Simulation 1 (Figure 15, 16 & 17). 

To check the simulations, wc also plot the V band 7?2i(=^2/^i) parameters (Simon & 
Lee 1981), resulting from Simulation 1, as function of log(P) in Figure 18. The filled circles 
are the original data and the crosses are the simulated data. Those crosses with triangles 
are the simulated data when either A\ or is about 2.0(7 away from the original values. 
It can be seen from the figure that the simulated data do indeed trace the original data. In 
addition, the P2/-P0 — 0.5 resonance around 10 days can be clearly seen from this figure. 

From Table 5 & 6, it can be seen that the errors of the means and each Fourier amplitudes 
are consistent in all three simulations. As expected. Simulation 2 gave the largest errors 
because the Gaussian noise generated in this simulation is bigger that other two simulations. 
Simulation 3 has the smallest errors, because this simulation is performed for one star. The 
errors for Simulation 1 is in between these two cases, which is typical. In addition, Table 
5 & 6 shows that the mean offsets of the simulated data and original data arc very small. 
This can be seen in Figure 15-17, as all histograms are almost centered at zero. Therefore, 
no bias introduced by the reconstruction procedure, because all four Fourier amplitudes are 
tightly clustered around the "true" values given by the original well sampled Cepheid. 

The simulations suggests that with typical HST Cepheid data, it is meaningful to fit a 
fourth order Fourier expansion and that our direct Fourier fitting procedures for estimating 
V and / band means are not biased. 

6. Discussion and Conclusion 

The light curve of a Cepheid can be reconstructed with the Fourier techniques described 
in this paper. For OGLE LMC Cepheid data, we have demonstrated that we can frequently 
improve the quality of the reconstructed light curve by restricting the range of Fourier 
amplitudes can take in the Fourier fit. These ranges are obtained from well sampled Cepheid 
light curves in the Galaxy and Magellanic Clouds. 
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We have then apphed these techniques to HST Cepheid data. Such data typically 
consist of 12 and 4 V and I band points respectively. In particular, Cepheids with sufficient 
observations permit a 4*^^ order Fourier expansion. However, the Cepheids in HST I band 
observations that have insufficient data points require the application of Fourier interrelations 
to reconstruct the 1 band light curves from the V band. In summary, the Fourier techniques 
for light curve reconstruction procedures, as applied to the HST data, include two main 
parts: 

• 4*'* order Fourier expansion to 12 V band data points. This is reconstructing the hght 
curves by a direct fit to the data points, with appropriate ranges for each of the Fourier 
parameters. 

• Fourier interrelation to 4 I band data points. This is using the linear relations of 
Fourier parameters between V and I bands to reconstruct the light curves. 

In general, these relatively simple techniques can reconstruct the Cepheid light curves 
quite well, given that the data points are relatively uniform and well sampled. These tech- 
niques serve as an alternative method to obtain the mean magnitude of Cepheids besides 
phase weighted intensity mean and template fitting procedures. However, the detailed com- 
parisons between these methods are beyond the scope of this paper. Some examples of the 
reconstructed light curves for extra-galactic Cepheids (observed by HOKP) with these Fourier 
techniques are presented in Figure 19, which show good agreements to the observed data 
points. 

We also performed Monte-Carlo simulations to determine the errors associated with 
our light curves reconstruction procedures. The errors of the mean magnitudes and Fourier 
amplitudes arc determined from the fits of Gaussian distribution to the error histograms. 
These Gaussian fits show that our procedure is unbiased. The errors are: 

• V band: amean = 0.034; a^i = 0.039; (Ta2 = 0.033; aAs = 0.028 and (Ta4 = 0.023. 

• I band: amean = 0.035; (Tai = 0.025; (7^2 = 0.021; ctas = 0.017 and (7^4 = 0.014. 

Could it be the case that the distribution of 12 points is so severe that a decent fit to 
the V band data is not possible? Certainly and it could be that in this case a template type 
technique will produce a V and I band mean. However, we again contend that in this case, 
it will probably be the case that this point source is not used in the final analysis. 
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6.1. The Correlations of Fourier Amplitudes 

It has been suggested by anonymous referee that the higher order Fourier amphtudes 
shown in Figure 3 & 5 are correlated with the first order Fourier amphtude (Ai) in the same 
bands (the Fourier intrarelations) . These have been examined, for example, by (Antonello 
et al. 1987) and in the appendix. However, the Fourier intrarelations are less preferable than 
Fourier interrelations for reconstructing the Cepheid light curves. The reasons are discussed 
in following examples and in Appendix. 

Here, we reexamine the results of Antonello et al. (1987) with our data sets and compare 
them to the Fourier interrelations. We used all "calibrating set" Cepheids and pick only the 
"good" OGLE LMC Cepheids (crosses) from Figure 3. Then we re-plot the Ai — A2 Fourier 
intrarelations and the Ai{V) — Ai{I) Fourier interrelation in Figure 20. The "Calibrating 
set" Cepheids and OGLE LMC Cepheids are in left and right panels, respectively. In the 
figure, we distinguished Cepheids with different period ranges: crosses are for Cepheids with 
periods shorter than 8 days (short-period Cepheid); triangles are for Cepheids with period 
in between 8 to 14 days (bump Cepheid); and filled circles are for Cepheids with periods 
longer than 14 days (long-period Cepheid). Error bars for each data point are omitted in 
the figure for clarity. The top two panels in Figure 20 are Fourier intrarelations in V and I 
band respectively, and the bottom panel is the Fourier interrelations of Ai(V) and Ai(I). 

It is clear from the figure that the scatter of Ai — A2 intrarelations are larger than the 
scatter of ^i(V^) —Ai{I) interrelations. Furthermore, the short-period, bump and long-period 
Cepheids populate different regions in the plots of intrarelations, as is also seen in Antonello 
et al. (1987). In contrast. The tightness of correlations in the Fourier interrelations is clear, 
and less dependent on period distribution. These two properties make Fourier interrelations 
more applicable in reconstructing the light curves than Fourier intrarelations. 

6.2. Reconstructing Light Curves for Short Period Cepheids 

The relatively large ranges of the Fourier amphtudes for Cepheids with periods less than 
10 days (or log(P) < 1.0) is one reason for using the Fourier techniques to reconstruct the 
Cepheid light curves. From Figure 3 to 5, it can be seen that the distribution of the Fourier 
amplitudes at periods longer than 10 days (or log(P) > 1.0) show certain trends, which can 
be used to construct templates light curves as a function of period. However, the Fourier 
amplitudes for Cepheids with period less than 10 days do not show any obvious trends but 
scatter around certain ranges (as given in Table 1). For example, at a given short period, 
Ai(V") may occupy the range from ~ 0.1 to ~ 0.4. Therefore, extra care has to be taken 
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when applying template fitting to the Cepheids with periods less than 10 days. In fact, at 
periods shorter than log(P) < 0.85 the template techniques do not apply (Stetson 1996) but 
the Fourier techniques still holds good. 

Although the extra-galactic Cepheids discovered in the past HST observations generally 
have period longer than 10 days, some short period period Cepheids have been discovered 
in Local Group galaxies IC 1613 (Dolphin et al. 2001) & Leo A (Dolphin et al. 2002) with 
HST and ground based observations, respectively. Hence for observing Cepheids at short 
period, the Fourier techniques are still useful to reconstruct the light curves. With the new 
installation of ACS (Advanced Camera for Survey) in HST, the discovery of more short 
period Cepheids in other galaxies is probable. 

6.3. The Effect of Metallicity 

One motivation for studying direct Fourier techniques to reconstruct Cepheid light 
curves is to deal with the possible metallicity dependence in Cepheids. Antoncllo et al. 
(2000) shown that the metallicity may affect the shapes of light curves close to 10 days, 
based on the comparison of Cepheid light curves in two galaxies. The appropriate ranges of 
the Fourier parameters (as given in Table 1) are assumed to cover the possible ranges due 
to the different metallicity environments. In addition, the Fourier interrelations are shown 
to be only weakly dependent on the metallicity though this does not mean that the actual 
light curve shape is independent of metallicity. 

The template light curves defined by Stetson (1996) are obtained from the Galaxy 
{Z = 0.020), LMC {Z = 0.008) and SMC {Z = 0.004), with sample of ~ 30-45 Cepheids in 
each galaxies. Stetson (1996) looked for and failed to find differences between the Galactic, 
LMC and SMC Cepheids in terms of their amplitude and light curve shape sufficient to 
change the estimation of mean magnitudes. However, much more data is available now. 
Paczyhski & Pindor (2000) found statistically significant different amplitudes between the 
Galactic/LMC/SMC Cepheids in the period range 1.1 < log(P) < 1.4, in the sense that 
higher metallicity Cepheids have higher amplitudes. Moreover, in the period range 0.0 < 
log(P) < 0.95, the amplitudes ratio of R21 and P31 for Cepheids in LMC (Andreasen & 
Petersen 1987) and SMC (Andreasen 1988) are found to be larger than the Galactic Cepheids. 
This may due to the different metaUicity environments, at least for the case of the SMC 
(Buchler & Moskalik 1994). Further some established workers in the field have found varying 
PC relations between the Galaxy, LMC and SMC and subsequently different PL relations 
(Tammann et al. 2001; Tammann & Reindl 2002). Thus while it may be that template 
methods yield a sufficiently accurate mean in varying metallicity environments, we contend 
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that the issue of light curve shape and metaUicity needs to be revisited in the hght of the 
data now currently available. 

6.4. Applications in Determining the Cepheid Distances 

An immediate application of these Fourier techniques is to reconstruct the Cepheid 
hght curves in NGC 4258 and derive its distance modulus. NGC 4258 is a spiral galaxy that 
has an accurate geometrical distance of 7.2 ± O.SMpc (corresponding to distance modulus 
of 29.28 ± 0.09) from water maser measurement (Herrnstein et al. 1999). By using the 
published data of 15 Ccphcids discovered from HST observations (Newman et al. 2001) 
and following the HOKP procedures (including using the same PL relations and reddening 
correction), wc derived a distance modulus of 29.38 ± 0.06 (random errors only) to NGC 
4258, which is about l.lcr away from water maser distance. If we apply the —O.OTmag 
correction due to the WFPC2 cahbration of "long vs short" exposures (this is caused by 
the charge-transfer efficiency of WFPC2, see HOKP and reference therein for details), then 
the final distance modulus becomes 29.31 ± 0.06, which is consistent with the water maser 
distance. In addition, Kanbur et al. (2002) have also derived the Cepheid distance to NGC 
4258 with the same techniques presented here but slightly different PL relations, the result 
of 29.36 ± 0.06 (random errors only) is also consistent to the water maser distance and the 
distance derived in this paper. 

When we apply these techniques to determine Cepheid distances to 16 galaxies using 
published photometry by the HOKP, our distance moduli, on average, are very similar to 
HOKP results, with a difference of ~ O.Olmag (Kanbur et al. 2002). On the other hands, 
application of these techniques to 3 Sandage-Tammann-Saha galaxies (NGC 3627: Saha et 
al. (1999); NGC 3982: Saha et al. (2001a); NGC 4527: Saha et al. (2001b)) yields shorter 
distance moduli of 0.07mags on average (Kanbur et al. 2002). The good agreements between 
the results in Kanbur et al. (2002) and in both HOKP and Sandage-Tammann-Saha galaxies 
suggest that our method is not too far of base. 

This work has been supported by HST Grant AR08752.01-A. Part of SN work was 
performed under the auspices of the U.S. Department of Energy, National Nuclear Security 
Administration by the University of California, Lawrence Livermore National Laboratory 
under contract W7405-Eng-48. We thank an anonymous referee for many helpful suggestions 
for making the paper more relevant. 
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A. 



The Fourier Intrarelations 



The technique of Fourier interrelations, as described in Section 4, are the hnear relations 
between the Fourier parameters in equation (1) in different wave bands. Similarly, there exist 
another group of linear relations between the Fourier parameters in the same wave band, 
known as the Fourier intrarelations. The reason for developing the Fourier intrarelations 
is same as in Fourier interrelation, mainly to reconstruct the I band light curves that only 
contain 4 epochs. Therefore, the Fourier intrarelations are the linear relation between the 1** 
order Fourier parameter and the higher order Fourier parameters in one particular wave band, 
since 4 data points only permit the 1** order Fourier expansion. The Fourier intrarelations 
have the following expression for either V or I bands: 



As in the case of Fourier interrelations, the coefficients are determined from the "cal- 
ibration set" Cepheids. The results of the fits to the data with equation (Al) & (A2) are 
presented in Table 7, which only fisted out the S^'' and 4*^* order Fourier intrarelations. The 
2nd Qj-(jgj- Fourier intrarelations were omitted because for sparse, 12 epoch V band data, 
we can fit the data with 2"*^ order Fourier expansion and expand to 4*'* order with Fourier 
intrarelations to reconstruct V band light curves. Then we can use the Fourier interrelations 
to reconstruct the I band light curves from the V band light curves. 

The plots for the Fourier intrarelations in "calibrating set" Cepheids are presented in 
Figure 21 & 22 for Fourier amplitudes and phases, respectively. From these figures, though a 
relation clearly exists, it may not be linear though we show the best fit linear relation. In the 
plots of A3 against Ai for both V and I band, there are some stars which lie well below the 
best fit linear relation. The A4 verses Ai plots also show some evidence of a non-linearity. 
Since the error bars have been plotted on these diagrams, the trends described here are real. 
This non-finearity may due to differences in long and short period Cepheids. 

As in the case of Fourier interrelations, the Fourier intrarelations in OGLE LMC 
Cepheids have also been found, and presented in Table 8. The corresponding plots of the 
Fourier amplitudes and phases are given in Figure 23 & 24, respectively. The possible non- 
finearity of the Fourier intrarelations that seem on the "calibrating set" Cepheids also show 
up in these figures, as some stars which lie below the best fit linear regression {A3 vs. Ai) 
and the indications of non-linearity in the A4 vs. Ai plots. Furthermore, by comparing 
the Fourier intrarelations in the "cafibrating set" and OGLE LMC Cepheids, we see clearly 



Aj = ttj + bjAi, j = 2, ...4 
(f)j = Cj + dj(pi, j = 2, ...4 



(Al) 
(A2) 
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that the slope of the best fit hne increases from the "calibrating set" Cepheids to the LMC 
Cepheids. The difference in slope between "calibrating set" and LMC Cepheids found in the 
intrarelations plots is real and is probably attributable to the metallicity differences. 

Due to the possible of non-linearity and metallicity dependence on the parent galaxy, 
the Fourier intrarelations are less preferable than the Fourier interrelations (Section 4) for 
reconstructing the I band light curves with only few data points available. 
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Table 1: The adopted ranges for the Fourier amphtudes, determined from Figure 3-5 for both 
V and I bands. 



Period Ranges Ai{V) A2{V) AsiV) A^{V) 



0.0 


< 


log(P) 


< 


1.0 


0.0- 


-0.40 


0.0- 


-0.18 


0.0- 


-0.11 


0.0- 


■0.08 


1.0 


< 


log(P) 


< 


1.1 


0.0- 


-0.44 


0.0- 


-0.14 


0.0- 


-0.08 


0.0- 


■0.07 


1.1 


< 


log(P) 


< 


1.2 


0.0- 


-0.48 


0.0- 


-0.16 


0.0- 


-0.10 


0.0- 


•0.08 


1.2 


< 


log(P) 


< 


1.3 


0.0- 


-0.52 


0.0- 


-0.19 


0.0- 


-0.13 


0.0- 


•0.09 


1.3 


< 


log(P) 


< 


1.4 


0.0- 


-0.52 


0.0- 


-0.23 


0.0- 


-0.16 


0.0- 


•0.10 


1.4 


< 


log(P) 


< 


1.5 


0.0- 


-0.52 


0.0- 


-0.26 


0.0- 


-0.16 


0.0- 


•0.10 


1.5 


< 


log(P) 


< 


1.6 


0.0- 


-0.49 


0.0- 


-0.25 


0.0- 


-0.14 


0.0- 


•0.10 


1.6 


< 


log(P) 


< 


1.7 


0.0- 


-0.48 


0.0- 


-0.23 


0.0- 


-0.13 


0.0- 


■0.09 


1.7 


< 


log(P) 


< 


1.8 


0.0- 


-0.43 


0.0- 


-0.20 


0.0- 


-0.12 


0.0- 


■0.08 


1.8 


< 


log(P) 


< 


1.9 


0.0- 


-0.40 


0.0- 


-0.17 


0.0- 


-0.11 


0.0- 


■0.07 


1.9 


< 


log(P) 


< 


2.0 


0.0- 


-0.35 


0.0- 


-0.13 


0.0- 


-0.10 


0.0- 


■0.06 



Period Ranges Ai{I) ^2(7) ^3(7) ^4(7) 



0.0 < 


log(P) 


< 


1.0 


0.0- 


-0.25 


0.0- 


-0. 


,12 


0.0- 


-0.08 


0.0-0.045 


1.0 < 


log(P) 


< 


1.1 


0.0- 


-0.28 


0.0- 


•0. 


,10 


0.0- 


-0.07 


0.0-0.04 


1.1 < 


log(P) 


< 


1.2 


0.0- 


-0.31 


0.0- 


•0. 


,11 


0.0- 


-0.08 


0.0-0.05 


1.2 < 


log(P) 


< 


1.3 


0.0- 


-0.33 


0.0- 


•0. 


,12 


0.0- 


-0.09 


0.0-0.06 


1.3 < 


log(P) 


< 


1.4 


0.0- 


-0.34 


0.0- 


-0. 


,13 


0.0- 


-0.10 


0.0-0.065 


1.4 < 


log(P) 


< 


1.5 


0.0- 


-0.34 


0.0- 


■0. 


,15 


0.0- 


■0.12 


0.0-0.07 


1.5 < 


log(P) 


< 


1.6 


0.0- 


-0.34 


0.0- 


-0. 


,16 


0.0- 


-0.11 


0.0-0.07 


1.6 < 


log(P) 


< 


1.7 


0.0- 


-0.31 


0.0- 


-0. 


,15 


0.0- 


-0.10 


0.0-0.065 


1.7 < 


log(P) 


< 


1.8 


0.0- 


-0.27 


0.0- 


-0. 


,13 


0.0- 


-0.09 


0.0-0.06 


1.8 < 


log(P) 


< 


1.9 


0.0- 


-0.24 


0.0- 


-0. 


,11 


0.0- 


-0.08 


0.0-0.055 


1.9 < 


log(P) 


< 


2.0 


0.0- 


-0.21 


0.0- 


-0. 


,10 


0.0- 


-0.07 


0.0-0.05 
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Tablc 2: The coefficients of the Fourier Interrelations, as determined from "calibrating set" , 
for fundamental mode Cepheids. Coefficients are corresponding to those defined in equation 
(3) & (4). Parameter characterizes the goodness of fit: = — a — Pxif/^ay^ + 



Relation a P "X? 

-0.006 ± 0.002 0.643 ± 0.007 2.69 x 10^ 

A2 0.001 ±0.001 0.600 ±0.011 9.74 xlO^ 

A3 0.000 ± 0.001 0.645 ± 0.023 6.05 x 10^ 

0.000 ± 0.001 0.631 ± 0.039 4.85 x 10^ 

Relation 7 ?7 

-0.178 ± 0.033 0.996 ± 0.001 2.21 x 10^ 

(/)2 -0.048 ±0.010 1.005 ±0.003 6.89 x 10^ 

03 -0.015 ±0.021 1.003 ±0.004 1.02 x 10^ 

04 -0.001 ±0.034 1.004 ±0.006 7.10 x 10^ 
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Table 3: The coefficients of the Fourier Interrelations for OGLE LMC (fundamental mode) 
Cepheids. The columns are the same as in Table 2. 



Relation a P 

A 0.004 ±0.005 0.601 ±0.002 2.64 x 10^ 

A2 0.001 ± 0.004 0.601 ± 0.004 1.68 x 10=^ 

^3 0.000 ± 0.001 0.606 ± 0.006 1.36 x 10^ 

^4 -0.001 ±0.001 0.625 ±0.015 1.16x10=^ 

Relation 7 ^ 

'Ji -0.208 ± 0.002 1.001 ± 0.001 4.79 x 10^ 

02 -0.067 ±0.004 1.001 ±0.001 1.49 x lO''^ 

03 -0.029 ±0.009 1.001 ±0.001 1.20 x 10^ 

04 -0.063 ±0.022 1.008 ±0.003 9.15 x 10^ 



Table 4: The coefficients of the Fourier Interrelations for OGLE SMC (fundamental mode) 
Cepheids. The columns are the same as in Table 2. 



Relation a P 

~Ai 0.003 ±0.001 0.608 ±0.002 3.87 x 10^ 

A2 -0.000 ±0.001 0.623 ±0.003 2.91 x 10^ 

^3 -0.001 ±0.001 0.650 ±0.005 2.65 x 10^ 

^4 -0.002 ±0.001 0.647 ±0.010 2.28 x 10=^ 

Relation 7 ^ 

'Ji -0.214 ±0.003 1.000 ±0.001 5.46 x 10^ 

02 -0.075 ± 0.004 1.001 ± 0.001 3.06 x 10^ 

03 -0.018 ±0.008 1.000 ±0.001 2.08 x 10^ 

04 -0.089 ±0.016 1.008 ±0.002 1.54 x 10^ 
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Table 5: The results of error analysis from simulations in V band. The Gaussian parameters 
II and (7 are mean off-sets and errors, respectively, for the means and the Fourier amplitudes. 



Simulation 


f^mean 




/^A2 




/^A4 


1 


0.0014 


-0.0121 


-0.0063 


0.0002 


0.0063 


2 


0.0006 


-0.0174 


-0.0062 


0.0013 


0.0095 


3 


-0.0018 


-0.0226 


-0.0154 


-0.0071 


0.0059 


Simulation 






crA2 




(^A4 


1 


0.0312 


0.0355 


0.0330 


0.0259 


0.0210 


2 


0.0337 


0.0393 


0.0332 


0.0280 


0.0228 


3 


0.0256 


0.0321 


0.0272 


0.0212 


0.0172 
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Table 6: The results of error analysis from simulations in I band. The Gaussian parameters 
II and (7 are mean off-sets and errors, respectively, for the means and the Fourier amplitudes. 



Simulation 


f^mean 




fJ'A2 


^A3 




1 


0.0012 


-0.0071 


-0.0053 


0.0018 


0.0048 


2 


0.0010 


-0.0109 


-0.0052 


0.0020 


0.0064 


3 


-0.0042 


-0.0132 


-0.0039 


0.0018 


-0.0004 


Simulation 






crA2 




(^A4 


1 


0.0283 


0.0226 


0.0195 


0.0158 


0.0130 


2 


0.0348 


0.0253 


0.0206 


0.0175 


0.0145 


3 


0.0251 


0.0208 


0.0163 


0.0156 


0.0108 
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Table 7: The coefficients of the Fourier Intrarelations, as determined from "cahbrating set" , 
for fundamental mode Cepheids. Coefficients are corresponding to those defined in equation 
(Al) & (A2). parameter is same as in Table 2. 



Relation 


a 


h 




A,{V) 


-0.043 ±0.002 


0.281 ±0.006 


3.97 X 10^ 




-0.021 ±0.002 


0.139 ±0.006 


2.97 X 10^ 




-0.025 ±0.002 


0.275 ±0.010 


2.77 X 10^ 


A,{I) 


-0.014 ±0.002 


0.145 ±0.010 


1.98 X 10^ 


Relation 


c 


d 




MV) 


2.544 ±0.020 


3.017 ±0.008 


2.42 X 10^ 


MV) 


0.202 ± 0.034 


4.061 ± 0.013 


1.40 X 10^ 


Ui) 


3.122 ±0.027 


3.043 ±0.011 


1.72 X 10=^ 


Ui) 


0.887 ± 0.044 


4.122 ±0.020 


9.02 X 10^ 
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Table 8: The coefficients of the Fourier Intrarelations for OGLE LMC (fundamental mode) 
Cepheids. The columns are same as in Table 4. 



Relation 


a 


b 




A,{V) 


-0.051 ±0.001 


0.364 ±0.005 


2.53 X 10^ 




-0.034 ± 0.002 


0.214 ±0.005 


1.71 X 10^ 




-0.031 ±0.001 


0.355 ± 0.002 


1.42 X 10^ 
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Fig. 1. — An example of the reconstructed light curve to a Cepheid in NGC 1326A. The 
solid curve is the reconstructed light curve from the least-square solution. It is clear that the 
least-square solution for this Cepheid does not provide a satisfactory fit to the data. Rather 
than solving for the parameters in equation (2), we fit the parameters in equation (2) with 
appropriate choice for the ranges (as given in Section 2.2) for Fourier amplitudes, Aj. The 
resulting light curve is drawn with a dashed curve, for comparison. The original data are 
indicated with error bars. 
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Fig. 2. — Two typical examples for the constructed light curves (in V band) by using the 
simulated annealing method for LMC Cepheids with similar periods. The ranges for the 
Fourier phases that are used in fitting the data points are from to 27r and the ranges for the 
Fourier amplitudes are: Ai : — 0.7; A2, A3, A4 : — 0.5. These choice of ranges are relatively 
large, as demonstrated in Section 2.1. Left: (a) Well constructed light curve. The data points 
are distributed relatively uniform in phase. This LMC Cepheid has log(P) = 0.4576. Right: 
(b) Poorly constructed light curve due to the bad phase coverage, and exhibits a numerical 
bump at phase ~ 0.5. This LMC Cepheid has log(P) = 0.4580. 



-31 - 



Fig. 3. — The distributions of the Fourier amphtudes for OGLE LMC Cepheids in V band. 
Total number of the Cepheids is 758, which includes 648 Cepheids with well constructed 
light curves (crosses) and 110 Cepheids with poorly constructed hght curves (triangles). 
The Cepheids in the "cahbrating set" are shown as filled circles, for comparison. 
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Fig. 4. — The distributions of the Fourier amphtudes for OGLE LMC Cepheids in V band 
(crosses) with well constructed light curves. The Cepheids in the "calibrating set" are shown 
as filled circles, for comparison. 
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Fig. 5. — The distributions of the Fourier amphtudes for OGLE LMC Cepheids in I band 
are shown as crosses. Some of the suspicious outUers are shown as triangles and labelled. 
The Cepheids in "calibrating set" are shown as filled circles, for comparison. 
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Fig. 6. — Examples of the constructed light curves for the outliers in Figure 5. The dashed 
lines are the light curves for Cepheid with similar periods, for comparison . 
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phase phase 

Fig. 7. — Examples of the reconstructed light curves fitted with the ranges of Fourier 
amplitudes given in Table 1. These light curves are indicated in solid curves. The 
dotted curves are the light curves fitted with initial ranges of Fourier amplitudes (i.e. 
Ao : - 0.7; ^1,^2, A3 : - 0.5). The periods of the OGLE LMC Cepheids are also 
given in upper-right corner. Original data points are also indicated. 
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Fig. 8. — Examples of the reconstructed light curves that exhibit numerical bumps/dips, 
regardless of fitting the data with smaller (solid curves, as given in Table 1) or larger (dotted 
curves) ranges of the Fourier amplitudes. This is mainly due to the bad phase coverage in 
the data, with gaps in between the data points. Original data points are also indicated. 
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Fig. 9. — The Interrelations among the Fourier amphtudes in the "cahbrating set" Cepheids. 
Crosses represent corresponding Fourier amphtudes and standard errors for each Cepheids. 
The dashed hues are the best-fit straight hues to the data, with the slopes and zero points 
given in Table 2. Solid and open squares are for bump and non-bump Cepheids. 
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Fig. 10. — Same as Figure 9, except for the Interrelations among the Fourier phases in the 
"cahbrating set" Cepheids. 
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Fig. 11. — The Interrelations among the Fourier amphtudes in the OGLE LMC (fundamental 
mode) Cepheids. The sohd lines are the best-fit straight hues to the data, with the slopes 
and zero points given in Table 3. The dashed lines are the best-fit model for "calibrating 
set" Cepheids, for comparison. 
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Fig. 12. — Same as Figure 11, except for the Interrelations among the Fourier phases in the 
OGLE LMC Cepheids. 
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Fig. 13. — The Interrelations among the Fourier amphtudes in the OGLE SMC (fundamental 
mode) Cepheids. The solid lines are the best-fit straight lines to the data. The dashed lines 
are the best-fit model for "calibrating set" Cepheids, for comparison. 
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Fig. 14. — Same as Figure 13, except for the Interrelations among the Fourier phases in the 
OGLE SMC Cepheids. 



Fig. 15. — Histograms of offsets for mean magnitudes in OGLE LMC Cepheids from light 
curve reconstruction procedures. The x-axes are the offsets of the simulated means from the 
original values, and the y-axes are the normalized counts. The dashed curves are the fitted 
Gaussian distribution to the histograms, with the parameters given in upper-left corners. 
Left and right panels are for V and I bands, respectively. 
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Fig. 16. — Histograms of offsets for OGLE LMC Ceplieids in light curve reconstruction pro- 
cedure in V band. Tlie daslied curves are tlie fitted Gaussian distribution to tlie liistograms, 
witli tlie parameters given in upper-left corners. 
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Fig. 17. — Histograms of offsets for OGLE LMC Ceplieids in light curve reconstruction pro- 
cedure in I band. Tlie dashed curves are the fitted Gaussian distribution to the histograms, 
with the parameters given in upper-left corners. 
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Fig. 18. — Plot of i?2i = ^2/^1 in V band for simulated data (crosses). The crosses with 
triangles are the simulated data points with either the values of Ai or A2 is ~ 2.0a away 
from the original values. The filled circles are original OGLE LMC Cepheid data. 




NGC 925, log(P)-1.515 




phase 



24 7 



NGC 7331, log{P)-1.463 



■ 24.5 r— f-— . 
25 r 
25.5 r 








phase 



23 L NGC 3351, log{P)-1. 633 
"23.5 F-- , ^^-^-"^ 



Fig. 19. — Examples of the reconstructed light curves in nearby galaxies. The V band light 
curves (left panel) are reconstructed via order Fourier expansion and the I band light 
curves (right panel) are reconstruct via the Fourier interrelations. 
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Fig. 20. — Comparisons of the correlations for Ai — A2 and ^i(V^) — Crosses are for 

Cepheid with period shorter than 8 days; triangles are for Cepheids with period in between 
8 to 14 days; and filled circles are for Cepheids with period longer than 14 days. Error bars 
are omitted for clarity. Left: "Cahbrating set" Cepheids; Right: OGLE LMC Cepheids. 
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Fig. 21. — The Intrarelations among the Fourier amphtudes in the "cahbrating set' 
Cepheids. The dashed hnes are the best-fit straight hnes to the data. 
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Fig. 22. — Same as Figure 21, except for the Intrarelations among the Fourier phases in the 
"cahbrating set" Cepheids. 
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Fig. 23. — The Intrarelations among the Fourier amphtudes in the OGLE LMC (fundamental 
mode) Cepheids. The sohd hues are the best-fit straight fines to tfie data. The dashed fines 
are the best-fit model for "calibrating set" Cepfieids, for comparison. 
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Fig. 24. — Same as Figure 23, except for the Intrarelations among the Fourier phases in the 
OGLE LMC Cepheids. 
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